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The progress in the modehng of exotic nuclei with an extreme neutron-to-proton 
ratio is discussed. Two topics are emphasized: (i) the quest for the universal mi- 
croscopic nuclear energy density functional and (ii) the progress in the continuum 
shell model. 



1. Introduction 

The goal of nuclear structure theory is to build a unified microscopic frame- 
work in which bulk nuclear properties (including masses, radii, and mo- 
ments, structure of nuclear matter), nuclear excitations (including a vari- 
ety of collective phenomena), and nuclear reactions can all be described. 
While this goal is extremely ambitious, it is no longer a dream. Indeed, 
hand in hand with experimental developments in the radioactive nuclear 
beam (RNB) experimentation, a qualitative change in theoretical modeling 
is taking place. Due to the infiux of new ideas and the progress in com- 
puter technologies and numerical algorithms, nuclear theorists have been 
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quite successful in solving various pieces of the nuclear puzzle. 

During recent years, we have witnessed substantial progress in many ar- 
eas of theoretical nuclear structure. The Effective Field Theory (EFT) has 
enabled us to construct high-quality NN and NNN bare interactions consis- 
tent with the chiral symmetry of QCD^'^. New effective interactions in the 
medium have been developed which, together with a powerful suite of ab- 
initio approaches, provide a quantitative description of light nuclei^'^'^'^'^. 
For heavy systems, global modern shell- model approaches''^' ^^'^^ and self- 
consistent mean- field methods^^'^^'^^ offer a level of accuracy typical of phe- 
nomenological approaches based on parameters locally fitted to the data. 
By exploring connections between models in various regions of the chart of 
the nuclides, nuclear theory aims to develop a comprehensive theory of the 
nucleus across the entire nuclear landscape. 

From a theoretical point of view, short-lived exotic nuclei far from sta- 
bility offer a unique test of those aspects of the many-body theory that 
depend on the isospin degrees of freedom^^. The challenge to microscopic 
theory is to develop methodologies to reliably calculate and understand the 
origins of unknown properties of new physical systems, physical systems 
with the same ingredients as familiar ones but with totally new and differ- 
ent properties. The hope is that after probing the limits of extreme isospin, 
we can later go back to the valley of stability and improve the description 
of normal nuclei. 

2. Towards the Universal Nuclear Energy Density 
Functional 

For medium- mass and heavy nuclei, a critical challenge is the quest for the 
universal energy density functional, which will be able to describe prop- 
erties of finite nuclei (static properties, collective states, large- amplitude 
collective motion) as well as extended asymmetric nucleonic matter (e.g., 
as found in neutron stars). Self-consistent methods based on the density 
functional theory (DFT) have already achieved a level of sophistication 
and precision which allows analyses of experimental data for a wide range 
of properties and for arbitrarily heavy nuclei. For instance, self- consistent 
Hartree-Fock (HF) and Hartree-Fock-Bogoliubov (HFB) models are now 
able to reproduce measured nuclear binding energies with an impressive 
rms error of ~700 keV 12,16,17 However, much work remains to be done. 
Developing a universal nuclear density functional will require a better un- 
derstanding of the density dependence, isospin effects, and pairing, as well 
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as an improved treatment of symmetry breaking effects and many-body 
correlations. 

2.1. Density Functional Theory and Skyrme HFB 

The density functional theory ^^'^^ has been an extremely successful ap- 
proach for the description of ground-state properties of bulk (metals, semi- 
conductors, and insulators) and complex (molecules, proteins, nanostruc- 
tures) materials. It has also been used with great success in nuclear 
physics^^'^^'^^'^^. The main idea of DFT is to describe an interacting sys- 
tem of fermions via its densities and not via its many-body wave function. 
The energy of the many body system can be written as a density functional, 
and the ground state energy is obtained through the variational procedure. 

The nuclear energy density functional appears naturally in the Skyrme- 
HFB theory^^'^^, or in the local density approximation (LDA)^^'^^, in which 
the functional depends only on local densities, and on local densities built 
from derivatives up to the second order. In practice, a number of local 
densities are introduced: nucleonic densities, kinetic densities, spin densi- 
ties, spin- kinetic densities, current densities, tensor-kinetic densities, and 
spin-current densities. If pairing correlations are considered, the number of 
local densities doubles since one has to consider both particle and pairing 
densities. 

In the case of the Skyrme effective interaction, as well as in the frame- 
work of the LDA, the energy functional is a three-dimensional spatial inte- 
gral of local energy density that is a real, scalar, time-even, and isoscalar 
function of local densities and their first and second derivatives. In the 
case of no proton-neutron mixing, the construction of the most general en- 
ergy density that is quadratic in one-body local densities can be found in 
Ref.^^. With the proton- neutron mixing included, the construction can be 
performed in an analogous manner^^. 

2.2. From finite nuclei to bulk nucleonic matter 

In the limit of the infinite nuclear matter, the density functional is reduced 
to the nuclear equation of state (EOS). The EOS plays a central role in 
nuclear structure and in heavy-ion collisions. It also determines the static 
and dynamical behavior of stars, especially in supernova explosions and in 
neutron star stability and evolution. Unfortunately, our knowledge of the 
EOS, especially at high densities and/or temperatures, is very poor. Many 
insights about the density dependence of the EOS, in particular the den- 
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sity dependence of the symmetry energy, can be obtained from microscopic 
calculations of neutron matter using realistic nucleon-nucleon forces^^'^^'^^. 
Those results will certainly be helpful when constraining realistic energy 
density functionals. Another constraint comes form measurements of neu- 
tron skin and radii^^'^^. Recently, a correlation between the neutron skin 
in heavy nuclei and the derivative of the neutron equation of state has been 
found^^'^^'^^, which provides a way of giving a stringent constraint on the 
EOS if the neutron radius of a heavy nucleus is measured with sufficient 
accuracy. 

A serious difficulty when extrapolating from finite nuclei to the extended 
nuclear matter is due to the diffused neutron surface in neutron-rich nuclei. 
As discussed in Ref.^^, the nuclear surface cannot simply be regarded as 
a layer of nuclear matter at low density. In this zone the gradient terms 
are as important in defining the energy relations as those depending on the 
local density. 



2.3. The First Step: Microscopic Mass Table 

Microscopic mass calculations require a simultaneous description of 
particle- hole, pairing, and continuum effects - the challenge that only very 
recently could be addressed by mean- field methods. A new development^^ 
is the solution of deformed HFB equations by using the local-scaling point 
transformation A representative example of deformed HFB calcula- 

tions, recently implemented using the parallel computational facilities at 
ORNL, is given in Fig. 1. By creating a simple load-balancing routine that 
allows one to scale the problem to 200 processors, it was possible to calcu- 
late the entire deformed even-even mass table in a single 24 wall-clock hour 
run (or approximately 4,800 processor hours). 

Future calculations will take into account a number of improvements, 
including (i) implementation of the exact particle number projection before 
variation^^; (ii) better modeling of the density dependence of the effective 
interaction by considering corrections beyond the mean-field and three- 
body effects^*^, the surface-peaked effective mass^^'^'', and better treatment 
of pairing^^; (iii) proper treatment of the time-odd fields and (iv) inclu- 
sion of dynamical zero-point fluctuations associated with the nuclear collec- 
tive motion ^3,44,45 ^^le density dependence is concerned, many 
insights can be obtained from the EFT^^. The resulting universal energy 
density functional will be fitted to nuclear masses, radii, giant vibrations, 
and other global nuclear characteristics. 
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Figure 1. Quadrupole deformations f3 (upper panel) and two-neutron separation ener- 
gies S2n in MeV (lower panel) of particle-bound even-even nuclei calculated within the 
HFB+THO method with Lipkin-Nogami correction followed by exact particle number 
projection. The Skyrme SLy4 interaction and volume contact pairing were used. (From 
Ref. 14.) 



Finally, let us remark that a realistic energy density functional does not 
have to be related to any given effective force. This creates a problem if 
a symmetry is spontaneously broken. While the projection can be carried 
out in a straightforward manner for energy functionals that are related to 
a two-body potential, the restoration of spontaneously broken symmetries 
of a general density functional poses a conceptional dilemma^^'^^. 



3. Continuum Shell-Model 

The major theoretical challenge in the microscopic description of nuclei, 
especially weakly bound ones, is the rigorous treatment of both the many- 
body correlations and the continuum of positive-energy states and decay 
channels. The importance of continuum for the description of resonances is 
obvious. Weakly bound states cannot be described within the closed quan- 
tum system formalism since there always appears a virtual scattering into 
the continuum phase space involving intermediate scattering states. The 
consistent treatment of continuum in multi-configuration mixing calcula- 
tions is the domain of the continuum shell model (CSM) (see Ref.^^ for a 
review). In the following, we briefly mention one recent development in the 
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area of the CSM, the so-caUed Gamow SheU Model. 
3.1. Gamow Shell Model 

Recently, the multiconfigurational CSM in the complete Berggren basis, the 
so-called Gamow Shell Model (GSM), has been formulated ^^'^^ The s.p. 
basis of GSM is given by the Berggren ensemble which contains Gamow 
states (or resonant states and the non-resonant continuum). The resonant 
states are the generalized eigenstates of the time-independent Schrodinger 
equation which are regular at the origin and satisfy purely outgoing bound- 
ary conditions. They correspond to the poles of the S matrix in the complex 
energy plane lying on or below the positive real axis. 

There exist several completeness relations involving resonant states 
In the heart of GSM is the Berggren completeness relation: 



where \un) are the Gamow states (both bound states and the decaying 
resonant states lying between the real A:- axis and the complex contour L^) 
and \uk) are the scattering states on L+. (For neutrons, / = resonances 
do not exist and, sometimes, one has to include the anti-bound / = state 
in the Berggren completeness relation s^'^^. This implies a modification of 
the complex contour which has to enclose this pole.) As a consequence 
of the analytical continuation, the resonant states are normalized according 
to the squared radial wave function and not to the modulus of the squared 
radial wave function. In practical applications, one has to discretize the 
integral in (1). Such a discretized Berggren relation is formally analogous 
to the standard completeness relation in a discrete basis of L^-functions and, 
in the same way, leads to the eigenvalue problem ^^|^) = E\^). However, as 
the formalism of Gamow states is non-hermitian, the matrix H is complex 
symmetric. 

One of the main challenges in the CSM is the determination of many- 
body resonances because of a huge number (continuum) of surrounding 
many-body scattering states. A practical solution to this problem has been 
proposed in Refs.^^'^^. It is based on the fact that resonances have signif- 
icant overlap with many-body states calculated in the pole approximation 
in which the Hamiltonian is diagonalized in a smaller basis consisting of s.p. 
resonant states only. The eigenstates representing the non-resonant back- 
ground tend to align along regular trajectories in the complex energy plane. 
As discussed in Refs. ^^'^^^ the shapes of these trajectories directly refiect 




(1) 
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the geometry of the contour in the complex /c-plane. In the two-particle 
case, this information can be directly used to identify the resonance states. 
However, this is no longer the case if more than two particles are involved. 

In the shell-model calculations with Gamow states, only radial ma- 
trix elements are treated differently as compared to the standard shell 
model. This means that the angular momentum and isospin algebra do 
not change in the GSM. However, expectation values of operators in the 
many-body GSM states have both real and imaginary parts. As discussed in 
Refs.^^'^^'^^, the imaginary part gives the uncertainty of the average value. 
It is also worth noting that, in most cases, the real part of the matrix ele- 
ment is influenced by the interference with the non-resonant background. 

Contrary to the traditional shell model, the effective interaction of GSM 
cannot be represented as a single matrix calculated for all nuclei in a given 
region. The GSM Hamiltonian contains a real effective two-body force 
expressed in terms of space, spin, and isospin coordinates. The matrix 
elements involving continuum states are strongly system-dependent and 
they fully take into account the spatial extension of s.p. wave functions. 

In the flrst applications of the GSM, a schematic zero-range surface 
delta force was taken as a residual interaction. As a typical example, the 
calculated level scheme of ^^O is displayed in Fig. 2 together with the 
selected E2 transition rates. It is seen that the electromagnetic transition 
rates involving unbound states are complex. 

The first applications of the GSM to the oxygen and helium isotopes look 
very promising^^'^^. The beginning stages of a broad research program has 
begun which involves applications of GSM to halo nuclei, particle- unstable 
nuclear states, reactions of astrophysical interest, and a variety of nuclear 
structure phenomena. The important step will be to develop effective finite- 
range interactions to be used in the GSM calculations. One would also like 
to optimize the path of integration representing the non-resonant contin- 
uum. In order to optimize the GSM configuration space, we intend to carry 
out GSM calculations in the Hartree-Fock basis. To this end, a Hartree- 
Fock program in the Gamow basis has been developed (GHF) The 
GHF method will also be applied to describe nuclear vibrational states in 
the continuum RPA (or QRPA) framework. 

4. Conclusions 

The main objective of this presentation was to discuss the opportunities in 
nuclear structure that have been enabled by studies of exotic nuclei with 
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Figure 2. The GSM level scheme of ■'^^O calculated in the full sd space of Gamow states 
and employing the discretized (10 points) non-resonant continuum. The dashed lines 
indicate experimental and calculated one-neutron emission thresholds. As the number of 
states becomes large above the one-neutron emission threshold, only selected resonances 
are shown. Selected E2 transitions are indicated by arrows and the calculated E2 rates 
(all in W.u.) are given (from Ref. ^^). 



extreme neutron-to-proton ratios. New-generation data will be crucial in 
pinning down a number of long-standing questions related to the effective 
Hamiltonian, nuclear collectivity, and properties of nuclear excitations. 

One of the major challenges is to develop the "universal" nuclear energy 
density functional that will describe properties of finite nuclei as well as 
extended asymmetric nucleonic matter as found in neutron stars. This 
quest is strongly driven by new data on nuclei far from stability, where new 
features, such as weak binding and altered interactions, make extrapolations 
of existing models very unreliable. 

Another major task is to tie nuclear structure directly to nuclear reac- 
tions within a coherent framework applicable throughout the nuclear land- 
scape. From the nuclear structure perspective, the continuum shell model 
is the tool of choice that will be able to describe new phenomena in dis- 
crete/continuum spectroscopy of exotic nuclei. 
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